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M. Olshanii [Phys. Rev. Lett. 81, 938 (1998)] recently solved the atom-atom scattering problem 
with a pseudopotential interaction in the presence of transverse harmonic confinement, i.e. within 
an 'atom waveguide', deriving an effective one-dimensional coupling constant that diverged at a 
"confinement induced resonance" (CIR). Here, we report numerical results for finite range potentials 
that corroborate this resonance. In addition, we now present a physical interpretation of this effect 
as a novel type of Feshbach resonance in which the transverse modes of the waveguide assume the 
roles of 'open' and 'closed' scattering channels. 
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Recently, there have been a number of experiments 
in which ultracold atoms and/or Bose-Einstein conden- 
sates have been loaded into magnetic or optical 'atom 

waveguides' [JES&SllSli. 0ne g° al of such 
experiments is to reach the 'single-mode' or quasi-lD 
regime, where only the ground state of transverse mo- 
tion is significantly populated at thermal equilibrium. 
This regime is of great practical interest due to the po- 
tential for ultra-sensitive rotation and gravitational gra- 
dient detection with guided single-mode atom interfer- 
ometers. In addition to such applications, reaching the 
quasi- ID regime is of significant theoretical interest as 
the ID delta-interacting boson gas represents one of the 
few known fully integrable quantum field theories. In 
a finite system with infinitely strong (hard-core) delta- 
function interactions, boson many-body states in ID have 
been shown to correspond via a one-to-one mapping with 
the highly-correlated states of the corresponding non- 
interacting Fermi gas . The properties of this Tonks- 
Girardeau gas have been a to pic of significant current 
theoretical interest |1 ES El El El Q El in antici- 
pation of future atom-waveguide experiments. Addition- 
ally, the homogeneous ID Bose with arbitrary-strength 
delta-function interactions, known as the Lieb Liniger 
model, is also a fully integrable system |T(| . 

To make the connection between experiments in tightly 
confining waveguides and theoretical models in ID, it is 
necessary to know the relationship between the effective 
ID coupling constant, gm, and the 3D scattering length, 
a. This problem was first addressed rigorously in |22| . 
where it was predicted that a 'confinement induced reso- 
nance' (CIR) modifies the effective interaction, resulting 
in an effective ID coupling strength which can be tuned 
from — oo to +oo by varying the transverse width of the 
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waveguide, a± over a small range in the vicinity of the 
resonance at a± = Ca, where C = 1.4603 . . .. Hence this 
resonance clearly has significant implications for atom- 
waveguide experiments. Until now, however, there has 
been no convincing physical explanation for the effect, 
thus raising questions concerning its appearance in sys- 
tems with finite-range interactions. 

The primary goal of this Letter, therefore, is to present 
numerical calculations of scattering in the presence of a 
cylindrical harmonic potential using finite-range atom- 
atom potentials, confirming the existence of the CIR. In 
addition, we provide a much-needed physical interpreta- 
tion of the effect as a Feshbach-type resonance involving 
bound states of the energetically closed transverse modes 
('channels'). Analogous numerical scattering studies in 
spherically symmetric traps have been reported |l7l Il8| , 
while scattering in the presence of harmonic confinement 
in one dimension is discussed in Isol j . 

We begin our analysis by considering a collision be- 
tween two atoms initially in the ground state of trans- 
verse motion. In the presence of harmonic transverse con- 
finement the center-of-mass motion and relative motions 
are separable, with the wavefunction of the relative coor- 
dinate satisfying an effective single-particle model with 
a stationary scatterer at the origin. Furthermore, if the 
longitudinal kinetic energy in the center-of-mass frame 
is less than the transverse level spacing then the atoms 
remain asymptotically frozen in the ground state. Low- 
energy scattering in this regime can then be modeled in 
the pseudopotential approximation [2lJ by the ID Hamil- 
tonian 

ft 2 d 2 

where z is the longitudinal atomic separation and p, is 
the reduced mass. In the pseudopotential approxima- 
tion, gm is typically obtained as in 

E3 by assuming 

that the wavefunction of the relative coordinate r = 
ri — T2 = zz + pp factorizes as VP(r) = (j>o(p)ip(z), where 
0o(p) = exp(— p 2 /2a 2 L )/(a±^yn) is the transverse ground 
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state, and a± = \Jfij pui±_ is the transverse harmonic os- 
cillator length for the relative atomic motion, lo± being 
the transverse trap frequency. The effective ID potential 
is then defined via / 27rpdp|0 o (p)| 2 ^^<5 3 (r) = g w S(z), 
where a is the 3D scattering length, which leads to 
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2h 2 a 
pa\ ' 



(2) 



From an exact solution of the three-dimensional scat- 
tering problem with a zero-range s-wave interaction, how- 
ever, it was recently predicted that a "confinement 
induced resonance" (CIR) modifies the effective interac- 
tion. By matching the low-energy scattering amplitude 
of the exact solution to that of Eq. it was found that: 
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where C = -((1/2) = 1.4603... We note that the use 
of Eqs. and © requires ka± <C 1, where hk is the 
collision momentum. The appearance of the resonance 
term in the denominator implies that to obtain an in- 
finite delta-function interaction (thereby accessing the 
Tonks- Girardeau regime), it is sufficient to satisfy the 
CIR condition a± ps Ca. From the naive formulation @ 
one would conclude that accessing this regime requires 
the more extreme condition a 3> m^, where n is the 
linear density. 

Before addressing the physical interpretation of the 
CIR, we first describe our numerical results pertaining to 
low-energy scattering and the effective ID coupling con- 
stant with finite-range interactions. The relation between 
giD and the ID scattering amplitude may be found by 
assuming that the scattering eigenstates ^(z) of Hamil- 
tonian take the form given in [2^ 



fe 



E = h 2 k 2 /2p 



(4) 



where f e is the coefficient for the even part of while 
the odd part can be shown to vanish by continuity ar- 
guments. From Eq. QJ, we find ~$>"(z) — — k 2 ^(z) + 
2ikf e (k)S(z), which must equal (2p,/h 2 )guj^(z) in the 
limit k — > 0. This gives 
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h 2 kf eR (k) 
k~b pf e i(k) 
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where f e R and / e j are the real and imaginary parts of 
f e . To extract analogous values for guj from numerical 
scattering calculations, we obtain eigenfunctions of the 
Hamiltonian H = H z + H± + V, where 



H -- h 92 ■ ft -- h 
z ~ 2pdz^ ±_ 2/i 



dp 1 



ld_ 

pdp 



+ ^V-(6) 



In the present work we restrict ourselves to the case 
of zero azimuthal angular momentum m=0, as there 
is negligible s-wave scattering for m ^ 0. For the 



atom-atom potential, we will study two cases: V(r) 
C 12 /r 12 -C 6 /r e (r 2 



z 2 + p 2 



and the spherical square- 
r), S(r) being the unit step func- 



well, V(r) = -VS(b 
tion. 

With application to the case of Cs atoms in a ID op- 
tical well |8j in mind, we consider atoms with the mass 
of 133 Cs, and with the Cq coefficients as determined re- 
cently 23] to be Cq (Cs) = 6890 a.u. For Cs atoms in the 
(F, M) = (3,3) state, the scattering length is determined 
only as an upper bound, a(Cs) < -140 nm, as compared 
with a_L = 29.5 nm reported in citeDWeiss. In order to 
study a more general situation, we allow C\i to vary. 
To simplify the numerics we consider the regime of just 
1 to 3 J = bound states, rather than the 47 bound 
states in the actual Cs2 3 S+ state, and we neglect other 
terms in the dispersion potential. The 6-12 potential may 
be characterized by R e , the minimum of the potential 
well, where R e = {2C 12 /C & ) 1 / 6 . For C 6 = C 6 (Cs) and 
for a series of values of Cn we obtain the free-space s- 
wave scattering length, a, by Numerov integration of the 
Schrodinger equation. When a resonance state passes 
through threshold a exhibits a simple pole, thus the 6-12 
potential provides the full range of 3D scattering lengths. 
For the second case of a spherical well potential the scat- 
tering length for a well of depth V and range b is given 
by a = b — t&n(br))/r), where rj 2 = 2pV/h 2 . 

To solve the scattering problem for such central poten- 
tials plus a transverse harmonic potential, we employed a 
numerical mesh in p and z, and found eigenfunctions in a 
cylindrical box of finite length. The box was sufficiently 
long in z that the asymptotic form of the wavefunction as 
\z\ — * oo could be determined. The eigenfunctions were 
of odd or even parity, with the odd parity functions ex- 
hibiting negligible scattering effects. For even functions 
of energy E = fijjj±_ + h 2 k 2 /2p : the asymptotic form is 



A/" [(1 + feR) cos (kz) - f eJ sin (k\z\)] Mp)- 

(7) 

Values for / 6 _r, f e i and AT were extracted from the coeffi- 
cients of sin (kz) and cos (kz) in Eq. Q, determined over 
a range of z values for which V(r) is negligible, and from 
conservation of probability current. From f e R and f e i, 
g±D is then obtained from Eq. (j3J) via extrapolation of 
finite k data to k = 0. The numerical mesh was provided 
by the discrete variable representation (DVR). In order 
to increase the density of points near z=0, the z coordi- 
nate was scaled by z = U(y), where U(y) — acosh(zfb), 
analogous to scaling used in [24|. A uniform mesh in y 
was then used along with a Laguerre DVR in p |26j| 
Eigenfunctions in first iteration were found by exact diag- 
onalization over a relatively small mesh (< 10,000 mesh 
points). As additional mesh points were added, sparse 
matrix diagonalization techniques were used. 

Numerical results for giu with the 6-12 potential and 
with the spherical well are shown in Fig. ^ in comparison 
with the analytic result (|3J). The numerical calculations 
clearly exhibit a singularity in gm at a/a± — 1/C . For a 
less than this value the calculated g\D values agree well 
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FIG. 1: The ID coupling constant, gm, in units of ft 2 /(/xox), 
as a function of a/aj_ for the 6-12 potential (stars) and the 
spherical well potential (triangles), as compared to the pseu- 
dopotential theory (solid line). The inset shows giD in units 
of h 2 /(n\a\) versus a±/a, to illustrate the behavior in the 
tight-confinement regime (a±/a — > 0). 



with the analytic expression. There do appear to be sys- 
tematic deviations between the exact numerical results 
and the pseudopotential theory that are greater for the 
6-12 potential than for the spherical well for a/a± > 1/C. 
This effect is connected to the Feshbach resonance in- 
terpretation and occurs because the bound state of the 
pseudopotential does not agree precisely with the up- 
per bound-state of the finite-range potentials having the 
same scattering length a. 

These numerical results validate the pseudopotential 
analytic result [22^. demonstrating that the CIR is in- 
deed a physical phenomena rather than an artifact of the 
pseudopotential approximation. It is therefore important 
to establish the physical origin of the CIR. In what fol- 
lows we are going to show that the CIR is in fact a zero- 
energy Feshbach resonance, occurring when the energy 
of a bound state of the asymptotically closed channels 
(i.e. the excited transverse modes) coincides with the 
continuum threshold of the open channel (lowest trans- 
verse mode). 

To illustrate this interpretation we consider the to- 
tal Hamiltonian H = H z + H± + V pseu do, where the 

pseudopotential operator is defined by (r| V pseu do\il)) = 
2-Kh 7 



^ -S 3 (r)-^r(r\tp) . Because our interpretation of the 
CIR relies heavily on the behavior of two-atom bound 
states, it is important to recall that the Huang-Fermi 
pseudopotential, V pse udo, supports a single bound state 
in free space with energy h 2 /(^,a 2 ), provided a > 0. With 
the existence of this bound state in mind, we proceed 
by formally splitting the Hamiltonian onto 'ground', 'ex- 
cited', and 'ground-excited coupling' parts according to 



H 



H„ + Hp 



H g -e 



where P n 



= PgHPg + P e HP e + (P e HP g +h.C.), (8) 



P e = X^^Li are the correspond- 



ing projection operators, |n) being the eigenstate of the 
transverse two-dimensional harmonic oscillator with ra- 



dial quantum number n and zero axial angular momen- 
tum (we recall that in denning H± we have already pro- 
jected out the states with non-zero axial angular momen- 
tum). The corresponding eigenvalues of the transverse 
Hamiltonian are given by H±\n, ) = hio±(2n + l)\n). 

The 'ground' Hamiltonian has a 1-d coordinate repre- 
sentation of the form of Eq. corresponding to the 
motion of a one-dimensional particle in presence of a 5- 
barrier with the 'bare' coupling constant given by Eq. 
(J2J . The spectrum of H g is continuous for energies above 
the threshold energy Ec, g = huj± (£ = 0) and contains a 
single bound state at energy Es, g < Ec, g for the case of 
positive 3D scattering length. Likewise, the spectrum of 
the 'excited' Hamiltonian is clearly continuous for ener- 
gies Ec, e = 3?kj^ {£ = 1), but as we will see below, H e 
supports one bound state of an energy Eb, s < Ec, e for 
all values of the 3D scattering length, a. 

It is clear that the remaining coupling between the 
'ground' and 'excited' subspaces will renormalizc the bare 
coupling J2| . Furthermore, according to the Feshbach 
scheme, one would predict a resonance in the renormal- 
ized gm for a set of parameters, such that the energy 
of the bound state of H e coincides with the continuum 



threshold of H„ 



Eb.c = En, 



CIR 



(9) 



As we will see below this scheme indeed predicts a posi- 
tion of the CIR exactly. 

The energy Eb, e of the bound state of H e can be found 
using the following two step procedure. First we iden- 
tify the bound state energy of the full Hamiltonian H 
as a pole of the scattering amplitude on the physical 
Riemann sheet. Second, we make use of the peculiar 
property of the two-dimensional harmonic oscillator that 
'excited' Hamiltonian H e &nd the full Hamiltonian H can 
be transformed to each other via a simple unitary trans- 
formation. This allows for simple relation between their 
energy spectra, in particular the energy of the bound 
states. 

The even-wave one-dimensional scattering amplitude 
f e at an energy E (E c a < E < E c , e ), as defined in 
has been derived in [22| . We have now obtained a closed- 
form analytic expression for the previously derived result, 
so that the scattering amplitude can be expressed as 



fe(k) 
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fca ± [^+C(l/2, -(ka ± /2)*)] 



(10) 



where 
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CO, a) = lim V 
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(JV + a) 



-s+l 



n^oo * — ' (n + a) 6 



S + 1 ' 

Re(s) > 0, Im{s) = 0, 

z s = | Z | S e ^(Arg( 2 )-2.) ; q < Arg(z) < 2?r (n) 



is the Hurwitz Zeta function 
k is given by E — Ec.g 
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and the wave-vector 
/c 2 /2/Lt. The bound state 
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FIG. 2: Numerical bound state energies and schematic of the 
Feshbach resonance scheme. The solid lines correspond to the 
analytic pseudopotential results. 

energies of the full Hamiltonian H will be given by the 
poles, fc, on the positive imaginary axis of the analytic 
continuation of f e (k): EbjuII = —h 2 Im 2 (k)/2fi. One 
can see that in order to avoid crossing the branch-cuts 
of the Zeta function, the continuation should be per- 
formed inside the < Arg(fc) < tt/2 quadrant of the 
complex plane. We find a single pole corresponding to 
the following implicit equation for the bound state en- 
ergy: C(l/2, -E Bifull /(2tkj ± ) + 1/2) = -o±/o. 

Notice now, that the full Hamiltonian H and the ex- 
cited Hamiltonian H e are connected via a simple trans- 
formation: H e = A 1 ' HA, where A' — J2^=o \ n + •"■>)( 7l l- 
Note that both H e and H include interactions, thus the 
above property is highly nontrivial and stems from the 
fact that the m = eigenfunctions of the two-dimensional 
harmonic oscillator all have the same value, l/y/lra±, at 
the origin. Thus the three dimensional ^-interaction has 
the same matrix elements between all the harmonic os- 
cillator states so that interaction matrix is unaffected by 
the shift operator. 

From the above we conclude that the bound state 
energy Es.e of the 'excited' Hamiltonian is related to 
the bound state energy Eb of the full Hamiltonian 
via Es.e = Eb + 2h~ujj_, and thus obeys an implicit 
equation ((1/2, -E B , e /(2hux) +3/2) = -a±/a. The 
CIR condition 10 can be now explicitly formulated as 
C(l/2, -E c , g /(2huj x ) +3/2) = -aja. Using E c , g = 
hu)j_i w e finally arrive at the exact CIR condition — = 
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C(l/2,1) = C(l/2,0) = -C. A similar effect is associ- 
ated with resonance behavior in harmonically confined 
2D scattering for a < 0|- We note that this reso- 
nance would most-likely be observed via changes in the 
macroscopic properties of the ground-state of a many- 
atom system, i.e. density distribution, as described in 

m 

This Feshbach scheme is illustrated in Figure 2, where 
we plot the bound-state energies Es, e (dark solid line) 
and Eb (thin solid line) as a function of the ratio a±/a. 
The continuum thresholds Ec, e an d Ec, g are also indi- 
cated, illustrating that the CIR occurs when the bound 
state of the manifold of closed channels, Es, e , crosses the 
continuum threshold of the open channel, Ec, g - In addi- 
tion, we have plotted the bound state energies of the full 
Hamiltonian as determined numerically for the 6 — 12 
and spherical well potentials, showing good agreement 
with the pseudopotential result. As the bound-state en- 
ergy deepens, we start to see quantitative disagreement 
between the bound state energies of the finite range po- 
tentials and the pseudopotential. This disagreement is 
consistent with the discrepancy in the position of the 
CIR shown in Fig. 1, showing that it is the bound state 
energy, and not the scattering length which determines 
the location of the CIR. Lastly, we note that while in free 
space a weakly bound state exists only for a > 0, we see 
that in the waveguide such a state exists for all a. These 
bound states may be of significant interest, allowing the 
formation of dimers via a modulation of the waveguide 
potential at the frequency (Ec,g —Es)/h. This may lead 
to an atom- waveguide based scheme for forming ultracold 
dimer molecules, as well as the possibility to use molec- 
ular spectroscopy as a sensitive probe of the atomic field 
inside the waveguide. 
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